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The three-dimensional, three-state Potts model is studied as a paradigm for high 
temperature quantum chromo dynamics. In a high statistics numerical simulation 



*vj ■ using a Swendson-Wang algorithm, we study cubic lattices of dimension as large as 

64 and measure correlation functions on long lattices of dimension 20 x 120 and 
j.^ ■ SO'^ X 120. These correlations are controlled by the spectrum of the transfer matrix. 

^H - 

t^ ■ This spectrum is studied in the vicinity of the phase transition. The analysis classifies 

^. 

^^ ■ the spectral levels according to an underlying 5*3 symmetry. Near the phase transition 

the spectrum agrees nicely with a simple four-component hamiltonian model. In the 

context of this model, we find that low temperature ordered-ordered interfaces nearly 

^ ■ always involve a disordered phase intermediate. We present a new spectral method 



for determining the surface tension between phases. 



12.38.Gc, 64.60. Cn, 12.38.Mh 



Typeset Using REVTEX 



I. INTRODUCTION 

The three-dimensional, three-state Potts model has long been studied as a paradigm 
for the phase structure of quantum chromodynamics (QCD) at high temperature in the 
heavy quark limit [U|. It has been found that at zero magnetic field, the Potts model has 
a weak first order phase transition, separating a low-temperature phase that breaks the 
Z(3) symmetry and a high-temperature phase in which the symmetry is restored 0. The 
underlying Z(3) symmetry requires that at low temperature there be three ordered or broken 
symmetry phases with the same free energy. In a finite volume system the behavior of the 
theory in the vicinity of the phase transition is complicated by tunneling among four phases: 
the three ordered phases and the disordered (symmetry restored) phase. 

The Potts model with no magnetic field corresponds to QCD with infinitely heavy 
quarks — in effect, without any dynamical quarks. A recent study of QCD without quarks 
by the APE group found that the correlation length of the system appears to grow as the 
physical volume is increased, suggesting an infinite correlation length in the infinite volume 
limit, a characteristic of a continuous phase transition |^. Doubts were raised that the 
phase transition is first order. The APE study differed from other contemporary work ^ in 
measuring correlation lengths on lattices with one long dimension. Two more recent high 
statistics studies with the Potts model and with QCD, using different indicators of the order 
of the phase transition, have reconfirmed the first order character of the phase transition in 
both models [^^. These more recent studies used a finite-size scaling analysis of susceptibil- 
ities to demonstrate a first order transition. The now widely suspected explanation for the 
confusion over correlation lengths is that finite volume tunneling among the Z(3)-equivalent 
ordered phases introduces a correlation related to the typical domain size of those phases. 
This correlation length does indeed become infinite in an infinite volume system, as a natural 
consequence of the first order character of the phase transition. 

To test this suspected explanation and to develop new insights into the effects of tun- 
neling, we have carried out a new study of the three-state, three-dimensional Potts model 



on lattices with one long j dimension [0. Our study emphasizes the determination of the 
spectrum of the transfer matrix. We demonstrate explicitly how tunneling modifies the 
spectrum and, as a consequence, the correlation lengths in the vicinity of the phase tran- 
sition. Our study parallels work done with the four- dimensional Ising model by Jansen et 
al [^ . Our work makes it possible to understand the results of the several references 0-|§] 
in a common framework. In the next section we discuss the S^ symmetry of the transfer 
matrix, develop a phenomenological four- component model for tunneling, and introduce the 
formulas needed for obtaining the spectrum. In Sec. |Tl|we present results of the simulation. 
We show that the spectrum and interface statistics agree well with predictions of the four- 
component model. We obtain the surface tension from the spectral splittings. In the final 
section we state our conclusions. 

II. PHENOMENOLOGY OF THE POTTS MODEL 

A. Transfer Matrix and Symmetries 

The three-dimensional, three-state Potts model is a classical spin system with one spin 
Si on each site i = {x, y, z) of a cubic lattice. Let the lattice dimension be L^ x Ly x L. 
Spins take on the values Si G Z{3) = {l,exp(±27ri/3)}. The partition function of the Potts 
model at zero magnetic field is given by 

Z{P) = Eexp l-PY.^es*s,+A , (1) 

Si V (iM) / 

where i + fi is the nearest neighbor site in the positive fi direction and the sum Sj is over all 
configurations of spins. Let the spins on a plane of constant lattice coordinate z be denoted 
Sz = {sx,y,z\x G 1, . . . ,Lx,y G 1, . . . ,Ly}. Then, as is well known, the partition function can 
be written as the trace of the transfer matrix raised to the power L: 

Z{(3)=TtT^, (2) 

where 



The transfer matrix can be regarded as a quantum mechanical operator acting on a Hilbert 
space of states described by configurations of spins arranged on the (x, y) lattice plane. All 
operators of present interest are local, and can be expressed as a function of the spins. The 
expectation value of an operator O on this ensemble is expressed in terms of the transfer 
matrix as 

(O) =Tr(T^O)/TrT^. (4) 

Local operators can be regarded as depending on a specific lattice plane z. The correlation 
between local operators A{z) and B{z') with z' > z is expressed as 

{A{z)B{z')) =Tt{T'AT''-'BT^-'')/TtT^. (5) 

A hamiltonian matrix H is defined in terms of the transfer matrix so that 

T = exp{-H). (6) 

Let us write the spectral decomposition of the transfer matrix in terms of the eigenvalues 
E„ of the hamiltonian H as 



■'n 



T = ^|n)exp(-K)(n|. (7) 

n 

Let n = denote the ground state. The correlation between two local operators A and B 
can be written in terms of the spectrum as follows: 

/ . . X Q. /x\ Emu {^\ A \n) {n\ B \m) exp[-(E„ - Em){z' - z)] exp{-EmL) 
<-'<'"'<' " - E„exp(-£„.) («' 

In this way correlations between local operators give information about the spectrum of 

the transfer matrix. (Of course, only the energy level differences {E^ — Eq) have physical 

significance.) 

The Potts model is symmetric under global transformations of the three-fold permutation 

group 5*3. These transformations are generated by 



5^^e±2W3g^ . s,^s* V,. (9) 

Therefore each of the eigenstates of H and each of the operators of interest O can be classified 
according to the three irreducible representations S (symmetric), A (antisymmetric) and M 
(mixed, two-dimensional). For example, the spin operator, itself, belongs to M, the operator 
|sp belongs to S, and the operator Re Si Im Sj — Im Si Re Sj belongs to A. 

B. A Simple Model 

To begin with the classification of states, consider first the extreme case f3 -^ oo. As is 
well known, the statistical ensemble reduces to three configurations of equal weight, with all 
spins aligned in either of the three Z(3) directions. In the Hamiltonian language the ground 
state of H at infinite (3 is three-fold degenerate, with all spins on the (x, y) plane aligned in 
one of the three directions. Call these three states |1), |2), and |3). These states are related 
by a Z(3) transformation R as follows: 

|2) = i?|l) |3) = i?|2). (10) 

At the other extreme [3 = the statistical ensemble contains all spin configurations with 
equal weight, and the ground state of the hamiltonian is not degenerate. Call it |0). These 
states have finite (3 counterparts. 

Beginning from these extremes, we introduce a phenomenological model for the spectrum 
at intermediate (3 in a finite volume system. As f3 is decreased from infinity, mixing between 
the three degenerate states occurs. Because of the S^ symmetry, the Hamiltonian matrix 
must be approximately of the form 



U e .^ 



e e 



v ' »/ 



(11) 



As a result of mixing, the degeneracy of the three states is lifted, giving rise to a symmetric 
ground state |0s) = |1) + |2) -|- |3) and a twofold degenerate mixed-symmetry state \0m) 



of slightly higher energy. The mixing parameter e depends on the transverse area L^Ly. 
Since mixing between the degenerate vacuums requires a rearrangement of the spins over 
the entire {x, y) plane, it is plausible that the dependence is 



exp{-pL.^Lya{p)), 



(12) 



where a{f3) is the surface tension for the interface. Indeed at large f3 it is easily shown that 
a{P) = 3/2. 

In the infinite volume system a phase transition takes place at /5 = jSc- In a finite volume, 
crossover occurs at j3c, where many observables change rapidly. For (3 < (3^ the ground state 
|0s) is then identified with the restored-symmetry-phase (disordered) vacuum |0), and the 
mixed symmetry state |0m) is identified as the lowest lying mixed-symmetry excitation of 
the symmetric vacuum. 

If the phase transition is first order, all four states coexist. Thus to model the rounding 
of the phase transition, we write the hamiltonian on the simplified basis of the four unmixed 
vacuum states, |0), |1), |2), and |3) as follows: 



' A A A A^ 



A e e 
A e e 
A e e 



(13) 



The row and column labels are in order 0, 1, 2, 3. The parameter A = a{(3 — f3c) with a > is 
the energy difference between the disordered (symmetric) phase state |0) and the degenerate 
ordered phase (broken symmetry) states |1), |2), |3). This difference vanishes at crossover. 
The parameter A gives the mixing strength between the disordered phase vacuum |0) and 
the degenerate ordered phase vacuums, and the parameter e gives the direct mixing strength 
between the ordered vacuums. The eigenenergies are 






+ -A - ivA2-4eA + 4e2 + 12A2 

2 2 



— e 



(14) 



El, = e + ^A + ^VA2-4eA + 4e2 + 12A2. 

As expected, there are two states belonging to the symmetric representation of S-^: one, the 
ground state |0s) and the other, an excited state |ls); and there is one twofold degenerate 
mixed symmetry state |0m) as before. For large /3 we expect the state |ls) to become degen- 
erate with a new excited mixed symmetry state |lr7i) (not included in our four- component 
model), in the same pattern as the |0s) and |0m) states. With our criteria for assigning 
transverse planes to the four phases (see Sec. P1B| below), the numerical simulation indi- 
cates that e << A. Thus in the region e << |A|, we may approximate the energy levels 
with 

1 . 1 



Eo. = 2^ - 2^^' + 12A2 (15) 

^om = (16) 

1 . 1 



Exs = 2^ + 2 v'A' + 12A2. (17) 

Figure |I| summarizes the expected behavior of the excitation energies, as a function of 
a(/? — /?c)/A, relative to the ground state energy, which has been renormalized to zero in 
this figure. The lowest two energy differences in this figure come from Eqs. ( |15| - |17|) and the 



energy difference E^^ is simply a sketch. We see that in a finite system the crossover results 
in a smooth connection between the energy levels on either side of /3c. In the infinite volume 
limit the crossover is much more rapid and the upper wings of the curves in Fig. |l| level off 
as a consequence of higher level crossings. 

By introducing a four-level system we have assumed a first order phase transition. In a 
continuous phase transition the states |0) and |0s) should be equivalent, and a three-level 
system would suffice. Although the model focuses on the lowest few levels, it can, of course, 
be enlarged to incorporate other excited levels as well. Indeed, in order to incorporate the 
excited levels quantitatively, at least eight levels are needed: four for the states already 
considered, and four more for the corresponding excited states. 

Let us be more explicit about the excitations of the ordered phase. In the infinite volume 



limit these vacuums are not mixed. We postulate that the lowest excitation of the state |1) is 
a state |1*) of mixed symmetry, reachable by acting on the vacuum with the zero momentum 
spin operator 

s(^) =IIs(^,j,,^). (18) 

x,y 

We shall often refer to the corresponding Schrodinger-picture operator s = s(0). Without 
loss of generality we introduce only one new state as follows: 

s|l) =7|1) + (5|1*). (19) 

The corresponding excitations in the other vacuums are reached by a Z(3) transformation 
R in analogy with Eq. ( PUJ ) as follows 

|2*) = i?|l*) |3*) = i?|2*). (20) 

Since RsR~^ = exp(— 27ri/3)s, the Z(3) symmetry requires 

s|2) =e2-'/3(7|2)+5|2*)) (21) 

s\3) =e-2'^*/3(7|3) + 5|3*)). (22) 

Now just as mixing between the degenerate vacuums leads to a symmetric state |0s) and two 
mixed symmetry states \Omk) {k = 1,2 labels the two degenerate components), we expect 
mixing among the excited state counterparts to result in a symmetric state |ls) and a mixed 
symmetry state llmk). With this notation we are implicitly identifying these states with 
the levels of Fig. |l[ From the 5*3 symmetry we then obtain explicit formulas for the matrix 
elements of the spin operator between these states: 



J2 |(ls|s|OmA;)|' =5^ 



k=l,2 

J2 |(0s|s|0mA;)|^ = 7^ 

fc=l,2 

Y, \{Os\s\lmk)f =6^ 

k=l,2 

Y: |(l.|.|lmA:)|^=7'. (23) 

k=l,2 
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We note also that if mixing among the degenerate states is weak the symmetric operator 
IsP satisfies 



'M |s|^ \0s) 



y^ (Im/cl |s|^ \Omk) 

k=l,2 



(24) 



C. Correlations 

We now write working formulas for the correlations between the spin operators s and 
|sp, based on the lowest lying levels discussed above. 

First, we observe that matrix elements must obey selection rules of the S3 symmetry. For 
example, the spin operator s is of mixed symmetry. The two components of the operator are 
just (Res,Ims). Thus matrix elements {n\ s \m) vanish if \m) and \n) are both symmetric. 
Since |sp is a symmetric operator, matrix elements {n\ |sp \m) vanish if \m) and \n) are not 
in the same 6*3 representation. 

Therefore, the leading nonvanishing terms in the correlation for |sp from Eq. (RI) are 



{\s{z)\^\s{0)\^) =\{0s\ \s\^ |Os)|' + |(l 



LS \S\ 



\0s) 



-Ei,,z 



+ e 



-Ei,{L-z)] 



{lm\ |sp |0m) 



'iEl„i~Eom)z 



+ e 



'iEl,n~Eoni)iL-z)T -EomL 



(25) 



and for s it is 



{s{z)s*{0)) = I (Om| s \0s) |2[e-^0'-^ + e~-^o,n(i-^)] ^ | ^^^| ^ |g^^ 



2\ -Ei^z I -Eim{L-z)-\ 



+ e 



+ I (Isl S \0m) |2rg-('E'is-So„)2 _j_ g-(£'i,-£;o„)(L-2)lg-£'omi 

+ I (lm| s \0m) |2[e-(i5i™-i?o™)2 ^ ^-(e^,^-Eo,^){l-z)i^^-Eo^l ^ 



(26) 
(27) 



The six transitions taken into account in these expressions are indicated in Fig. |^. 

Clearly, in order for the several transitions for be discernible in the fits to correlation 
functions, it is necessary that the spectral components be both strong and well-separated. 
With our data we are able to distinguish two spectral components for the mixed opera- 
tor and find only one significant spectral component for the symmetric operator. These 
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transitions correspond to dropping all terms with the factor exp(— E'om-^)- For the mixed 
operator, this approximation can be justified as follows: (1) Over the range (3 < (3c, the 
factor exp(— ii^om-^) is small (at the largest, of order 1/10). This factor multiplies the third 
term on the rhs in both expressions and the fourth term in the second expression. Appealing 
to our phenomenological model and Eqs. ( P^D and (P^, which makes it possible to compare 
the second and third terms, we see that we may drop terms in exp(— E'om-^) from the sum- 
mation, thereby eliminating all transitions leading to the level -Eom- (2) Over the range 
j3 > (3c, the spectral lines for all transitions from the levels Is and Im to Os and Om are too 
close to be resolvable with our statistics. Thus the second spectral component in the mixed 
operator is presumably a composite of all three transitions permitted by the selection rules. 
For the symmetric operator our failure to locate the second spectral component presumably 
refiects an insufficiently strong coupling to this operator. 

We are left, finally, with a three-parameter expression for the symmetric operator and a 
four-parameter expression for the mixed operator: 



{s{z)s*{0)) = Ciie-^o™" + e-^"'"^^-^)] + Caie"^!'"^ + e-^i-^^"^)]. (28) 



III. NUMERICAL SIMULATION 

Since the simulation of tunneling effects requires an algorithm that gives efficient and 
rapid sampling of the phase space, particularly in the crossover region, we used the 
Swendson-Wang (SW) algorithm [^. The code was checked by comparing measured observ- 
ables in selected extensive runs with results from simulations using a heat bath algorithm. 

Simulations were carried out on cubic lattices of size A^'^, for N = 16, 20, 24, 30, 48, 
and 64 in order to confirm the first order nature of the phase transition through the use of 
finite-size scaling analysis. For this purpose we measured the specific heat and the fourth- 
order cumulant. Table | summarizes the extent of this data sample. The majority of the 

10 



simulations were carried out on "cylindrical" lattices to obtain the spectrum of the transfer 
matrix. These lattices were of size 120 x 20^ for 18 values of (3 and with size 120 x 30^ for 
17 values of (3, as summarized in Table |T[ Runs of as many as two million sweeps made it 
possible to gain control of the correlated data and to obtain good estimates of the parameter 



errors. 



A. Simulations on cubic lattices 

In Fig. ^ is shown a histogram in the average energy for the 64^ lattice near (3^- A clean 
separation of the phases is apparent. The importance of finite size scaling for determining the 
order of a phase transition has been repeatedly emphasized |0,^. Although these references 
already provide excellent confirmation of the first order character of the phase transition, 
our high statistics results make an even stronger case. Indeed, with such a clean separation 
of phases in the 64'^ lattice, it is scarcely necessary to belabor the point. Nevertheless, we 
present the finite size scaling results for the sake of completeness. The fourth order cumulant 

W,^) = l-^, (29) 

is minimum near the crossover (3^- At a continuous phase transition, the minimum value 
Vl^Pc, V) tends to | in the infinite volume limit. At a first order phase transition, however, 
the limit is not so constrained, but obeys a scaling law, 

V,{VU„. - 1 - f^l^^'ll + OiyV), (30) 

where £"+ and E^ are the most probable energies in the two coexisting phases. 

The minimum value of the cumulant was found by combining measurements at a range 
of (3 values near /3c, using a Ferrenberg-Swendson "scanning" or "histogram" technique [ITT 



The resulting minimum values are plotted in Fig. |^. A linear fit yields the asymptotic value 
0.647(3), clearly distinct from |. 
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We turn now to the specific heat. At finite volume the specific heat C„ peaks at the 
crossover. The maximum value C^^max increases with increasing volume. If the phase tran- 
sition is first order, the peak scales as 

Cv^max = a + bV. (31) 

The maximum is again determined using the Ferrenberg-Swendson technique. These values 
are plotted in Fig. |^. A linear increase is apparent. 

A value of f3c can be fixed, either from the peak in the specific heat, or from the minimum 
of the cumulant Vl- The two values do not necessarily agree at finite volume, but should 
agree in the infinite volume limit. In the infinite volume limit we find (3c,oo = 0.36704(2) in 
agreement with Gavai, Karsch and Petersson 0]. 

B. Simulations on cylindrical lattices 

Measurements on the asymmetrical lattices were taken every 250 SW sweeps. Observ- 
ables recorded were these: the spin averages s{z) as a function of z [Eq. (^8|)], the average 
energy, 

E = PY.^es*s,+^/V, (32) 

and the number of clusters N^. Subsequent analysis produced the symmetric and mixed 
operator correlations and spectrum, the mean spins, and the projected-spin order parameter 

Sproj = max Re(s, e^"*/^s, e'^'^/^s), (33) 

where s = Y. Si/V. Also constructed were interface statistics. They are described below. 
For observables not discussed here, see Ref. [0]. 

a. Spectrum of the Transfer Matrix Correlations in the operators |s(2;)p and s{z) were 
measured and fit to the formulas (^) for both transverse sizes 20^ and 30^. As usual, fiuc- 
tuations in the measurements were strongly correlated in z, so it was necessary to determine 
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these correlations and incorporate them in the x^ analysis. Because of the large size of the 
data sample, it was possible to use all principal factors in the analysis of covariance. The 
spectrum was determined from a global fit to the data. The fitting range began at a mini- 
mum distance Zmm and extended to the full length of the lattice. The minimum distance was 
varied until a semblance of a plateau in the spectrum was reached, within the determined 
statistical errors. The values quoted are based on the minimum distance that gave the 
highest confidence level for the fit. The minimum distance thus determined varied smoothly 
from 2 for the mixed operator at the smallest (3 where the correlation length is shortest to 8 
at the largest (3 where the correlation length is longest. For the symmetric operator Zmin was 
5 for the smallest and largest /?'s, ranging gradually to 15 — 20 for the intermediate values 
where the correlation length is longest. The resulting spectrum is summarized in Figs. |^, |^, 
1^, and I and in Tables |T| and |^ There is obviously a strong resemblance with features of 



the simple model of Sec. |IIB| . From Fig. ^ and P, we see that the correlation length tends to 



infinity as (3 increases above /5c, just as with lattices of similar geometry in SU(3) Yang-Mills 
theory |^. As we have seen in the simple model, this feature is an expected consequence of 
finite volume tunneling between the degenerate ordered phases. 

b. Tunneling Statistics In Fig. |^ we plot the slice spin averages s{z) for a representative 
configuration at /5 > /?c, showing tunneling between the ordered phases. We devised two 
statistics: N^io, to give a measure of the number of phase boundaries between the disordered 
and one of the ordered domains, and a statistic Ngo, for phase boundaries between two 
ordered domains. A portion of the lattice was considered to be in the disordered phase 
on the plane z, if \s{z)\ < 0.23 for at least three consecutive values of z. If \s{z)\ > 0.23 
for at least three consecutive values of z, the lattice plane was considered to be in one of 
the three ordered phases, according to the value of args(z). The value 0.23 was chosen to 
correspond to the minimum of the histogram of occurrences of values of |s| at the crossover, 
and so corresponds to a value intermediate between the disordered and ordered phases. The 
requirement of three consecutive planes was adopted to permit occasional excursions from 
the ideal value of 0.23 within a single phase. We encountered no configurations that did not 

13 



have at least three consecutive planes. Indeed the phase coherence was extremely high with 
many lattices consisting of a single phase. 

Obviously our classification criteria are arbitrary. Our approach differs from that of 



Karsch and Patkos [|T2[ who classified all boundaries as type oo for /5 > (3c- The ambiguity all 
methods must deal with is distinguishing a broad interface between two ordered phases from 
a transition to an intermediate disordered phase. Any definition must recognize, however, 
that in perturbation theory, finite volume mixing necessarily produces a disordered phase 
intermediate for /5 slightly above Pc- 



Shown in Figs. |T0|, [11], and |1^ are results for the measure of the mean numbers Noo and 
Nfio- It is apparent that with our definitions of these boundaries, a transition between two 
ordered phases is unlikely to take place directly, but proceeds through what we identify as 
a disordered phase intermediate. (Figure |10| shows the number of oo phase boundaries for 
the 120 X 20^ lattices. The corresponding number for 120 x 30^ is negligible.) Thus if we 



interpret these results in terms of the model of Sec. [11 B|, we find that e << A. Thus the 



approximation, Eqs. (|T3|-|T^) applies. Inverting these expressions, the parameters A and A 
of the simple model can then be derived from the observed spectrum: 

A = ^1, - 2Eom (34) 



A = yEorn{Els — -Eom)/3 (35) 

Let us estimate the number of "do" (ordered/disordered) phase boundaries expected in 
our four-component model. Since A in the four-component hamiltonian mixes the disordered 
and ordered states, we can introduce a chemical potential for do boundaries by replacing A 
with Xe^. Thus the number of do boundaries is just 

(iVdo) = A^ In Z, (36) 

but from Eqs. (0), (|]), and (p!5|) we have InZ ^ — L£'os(A), so 



{Ndo) ~ 6LAVVA2 + 12A2. (37) 



14 



This predicted value (for (3 > /3c) is plotted together with the observed number in Fig.[TT. 
The agreement is quite satisfactory. 

c. Surface Tension The surface tension between two ordered phases can be estimated 
in perturbation theory (dilute interface approximation) through 

AVA = e-f^""^"^. (38) 

The quantity on the left is the transition probability between two ordered phases via a disor- 
dered phase intermediate in second order perturbation theory, based on the four- component 
hamiltonian, and the quantity on the right is the Boltzmann weight for the interface. Thus 
estimated, it is plotted in Fig. O. There is a strong dependence on the transverse size. 



suggesting a significant finite size correction. Karsch and Patkos [|T2|] measured the surface 
tension in the Potts model and also find a strong dependence on transverse size. They follow 
a statistical approach based on the frequency of oo interfaces, using a different definition 
from ours, which excludes all do interfaces. For example the configuration of Fig. ^ yields 
nine "bare domain walls" in their language, but all presumably classified in their scheme 
as "defects" and not genuine interfaces. Perhaps it is not surprising, therefore, that our 
surface tension is two to three times lower than theirs (after making allowances for differ- 
ent conventions in the hamiltonian). Clearly the determination of surface tension is model 
dependent. 

IV. SUMMARY AND DISCUSSION 

Our high statistics study of the three-state three-dimensional Potts model using the 
Swendson-Wang updating scheme once again confirms the first order character of the phase 
transition. At volumes as large as 64'^ the separation of coexisting phases is so clear that a 
sophisticated finite size scaling analysis is scarcely necessary. 

Exploiting an S3 symmetry in lattices with one long dimension, we have obtained the 
lowest spectral levels of the transfer matrix in this model and find excellent agreement with 
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the spectrum of a simple four-component model, featuring a first order phase transition. 
This analysis provides a clear explanation of the mechanism that gives rise to an infinite 
correlation length in the low temperature phase. 

The statistics of phase boundaries at low temperature are consistent with a perturbative 
treatment of tunneling in the four-component model. With our assignment of planes to 
phases, we find that ordered-ordered phase boundaries almost always involve an intermediate 
disordered phase. We introduced a spectral method for estimating the surface tension, based 
on the four-component model. 

One important goal of finite size spectral analysis is to remove tunneling-related finite- 
size effects from the spectrum, with the hope of extracting the infinite volume values of 
the excitation spectrum. Would this be feasible using our methods? Unfortunately, to 
remove significant finite size effects apparently requires introducing more parameters into 
the four-component model and into our fitting functions than data of the quality of ours 
warrants. The four component model would have to be augmented by at least four more 
components. Thus we must rely upon alternate, empirical methods. For example, Fukugita 
et al 1^ calculate in a cubic volume. In the vicinity of the phase transition, they classify 
configurations into two groups: disordered and ordered, according to the value of the global 
order parameter. They then measure "pure-phase" correlation lengths in each subset. For 
small volumes there is a region of overlap in which this classification risks misidentification of 
the phase. The contamination of incorrectly classified configurations decreases as the volume 
is increased and the overlap decreases. Thus one may hope for an empirically determined 
extrapolation to the infinite volume limit. 
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FIGURES 
FIG. 1. Phenomenological model of the lowest three energy levels at crossover. The ground 

state energy has been adjusted to zero at all (i. The level Eim does not come from the model; it is 



based merely on a guess. The others come from the four-component model Eq. (15|-17D. Transitions 



induced by a mixed operator are indicated with "M" and a symmetric operator with "S" . 

FIG. 2. Histogram in the average energy at /3 = 0.36705 (near the critical value) on a 64^ 
lattice, showing cleanly separated coexisting phases. 

FIG. 3. Finite size scaling of the minimum in the fourth order cumulant. 

FIG. 4. Finite size scaling of the peak in the specific heat. 

FIG. 5. Energy level Eis of symmetry S for the 20^ x 120 lattice, as a function of (3. 

FIG. 6. The same, but for the 30^ x 120 lattice 

FIG. 7. Energy levels i?om and Eim of symmetry M for the 20^ x 120 lattice, as a function of 

(5. 

FIG. 8. The same, but for the 30^ x 120 lattice 
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FIG. 9. Plot of the complex modulus and argument of the spin vector s(z) vs z (averaged over 
the transverse plane) for a typical configuration selected from the data sample at /3 = 0.3672 (near 
the phase transition) on a 20^ x 120 lattice. The plot symbol indicates the phase to which the lattice 
plane is assigned, based on modulus and argument. The vertical bars indicate an assigned phase 
boundary, based on our arbitrary rule that at least three consecutive planes must be classified in 
that phase. The horizontal line in the modulus plot indicates our division between ordered and 
disordered phases. Two ordered-ordered phase boundaries appear in this configuration, each with 
a disordered phase intermediate. 

FIG. 10. Average number of phase boundaries separating ordered symmetry phases for 
120 X 20^. 

FIG. 11. Average number of phase boundaries separating a disordered phase from an ordered 
phase for 120 x 20^. Crosses are measured directly. Octagons are calculated from the spectrum in 
the four-component model. 

FIG. 12. Same, but for 120 x 20^. 

FIG. 13. Ordered-ordered interface surface tension determined by spectral methods in the 
four-component model. Octagons are for 20^; crosses are for 30^. 
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TABLES 



TABLE L Data sample for cubic lattices 



V 



P 



sweeps (10^ 



16^ 
643 



0.3663 
0.3665 
0.3667 
0.3669 
0.3670 
0.36705 



0.6 
0.5 
0.7 
0.7 
0.5 
0.5 
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TABLE II. Data sample for cylindrical lattices. (Million sweeps.) 
p 120 X 20^ 120 X 30^ 

.3650 1.0 1.0 

.3655 — 1.0 

.3660 1.0 1.0 

.3665 1.0 1.0 

.3668 2.0 1.0 

.3669 2.0 2.0 

.36695 2.0 2.0 

.3670 2.0 2.0 

.36705 2.0 2.0 

.3671 2.0 2.0 

.36715 2.0 2.0 

.3672 3.0 2.0 

.36725 2.0 2.0 

.3673 1.0 1.0 

.3674 1.0 1.0 

.3675 1.0 1.0 

.3678 1.0 — 

.3680 1.0 1.0 

.3685 1.0 — 
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TABLE III. Effective masses and couplings for the 120 x 20^ lattice 



P 



A, 



A. 



E 



Is 



Ci 



E, 



Om 



Co 



Elr 



.3650 


2.38(1) 


0.5(1) 


.231(27) 


2.9(1) 


.167(5) 


0.18(11) 


.54(28) 


.3660 


3.20(3) 


0.94(7) 


.168(10) 


3.3(9) 


.10(2) 


1.0(9) 


.20(7) 


.3665 


4.36(4) 


1.6(2) 


.118(9) 


5.8(1) 


.072(1) 


0.31(4) 


.54(11) 


.3668 


6.11(6) 


2.6(3) 


.109(8) 


8.0(1) 


.048(1) 


0.30(5) 


.34(11) 


.3669 


7.17(7) 


2.5(2) 


.088(6) 


9.3(1) 


.0397(6) 


0.34(9) 


.41(12) 


.36695 


7.8(1) 


2.9(4) 


.095(11) 


10.1(2) 


.037(1) 


0.34(1) 


.41(18) 


.3670 


8.51(8) 


2.7(4) 


.088(10) 


10.8(1) 


.0346(6) 


0.24(6) 


.31(16) 


.36705 


9.59(9) 


2.8(4) 


.084(9) 


11.8(1) 


.0302(4) 


0.26(8) 


.38(17) 


.3671 


10.5(1) 


3.0(4) 


.084(9) 


12.4(1) 


.0272(4) 


0.35(5) 


.28(8) 


.36715 


11.8(1) 


2.8(4) 


.082(9) 


13.2(3) 


.0233(6) 


0.6(3) 


.12(5) 


.3672 


13.2(1) 


2.8(1) 


.082(3) 


14.1(4) 


.0207(6) 


0.5(3) 


.11(7) 


.36725 


14.4(1) 


2.6(2) 


.081(5) 


14.5(6) 


.0180(7) 


0.8(4) 


.10(5) 


.3673 


15.8(2) 


2.8(4) 


.085(11) 


15.4(4) 


.0159(7) 


0.8(2) 


.14(9) 


.3674 


18.6(2) 


1.8(2) 


.068(12) 


16.2(3) 


.0122(5) 


1.8(8) 


.21(7) 


.3675 


20.7(2) 


1.5(3) 


.076(14) 


16.9(3) 


.0102(5) 


0.9(3) 


.18(8) 


.3678 


25.5(1) 


1.3(1) 


.142(11) 


15.6(5) 


.0039(6) 


1.3(1) 


.14(3) 


.3680 


27.2(1) 


1.12(5) 


.175(9) 


15.9(4) 


.0029(4) 


1.2(2) 


.23(4) 


.3685 


30.4(1) 


1.03(6) 


.256(13) 


15.7(3) 


.0007(3) 


1.2(1) 


.28(3) 
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TABLE IV. Effective masses and couplings for the 120 x 30^ lattice 



P 



A, 



A. 



Els 



Ci 



E, 



Om 



Co 



E 



Im 



.3650 


1.02(1) 


0.2(1) 


.30(5) 


1.2(1) 


.182(8) 


.1(1) 


.39(16) 


.3655 


1.11(1) 


0.23(6) 


.25(3) 


1.29(8) 


.157(6) 


.15(7) 


.45(14) 


.3660 


1.23(1) 


0.38(8) 


.25(3) 


1.56(2) 


.143(3) 


— 


— 


.3665 


1.48(1) 


0.44(9) 


.15(2) 


1.74(9) 


.101(4) 


.24(9) 


.30(6) 


.3668 


2.03(3) 


0.74(6) 


.065(7; 


2.3(2) 


.052(3) 


.59(18) 


.16(3) 


.3669 


3.21(6) 


1.28(5) 


.061(4; 


3.7(1) 


.031(1) 


.72(12) 


.13(1) 


.36695 


4.09(9) 


1.61(4) 


.044(4; 


5.2(1) 


.0262(6; 


.39(4) 


.19(2) 


.3670 


6.3(2) 


2.02(7) 


.038(5; 


7.3(1) 


.0185(3; 


.41(3) 


.20(2) 


.36705 


9.9(3) 


2.1(1) 


.033(6; 


9.9(1) 


.0133(2; 


.38(3) 


.21(3) 


.3671 


14.4(2) 


1.6(1) 


.039(6; 


11.8(1) 


.0092(2; 


.55(4) 


.14(1) 


.36715 


17.4(3) 


1.3(2) 


.030(9; 


12.6(1) 


.0064(2; 


.65(3) 


.16(1) 


.3672 


19.6(1) 


0.98(3) 


.066(4; 


12.5(2) 


.0044(2; 


.75(4) 


.12(1) 


.36725 


20.9(1) 


0.85(8) 


.085(7; 


12.2(2) 


.0028(2; 


.86(3) 


.12(1) 


.3673 


21.6(1) 


0.7(1) 


.09(1) 


12.6(2) 


.0027(3; 


.73(4) 


.13(1) 


.3674 


23.1(1) 


0.5(2) 


.11(2) 


12.3(3) 


.0011(4; 


.69(4) 


.15(2) 


.3675 


24.1(1) 


0.5(1) 


.13(1) 


12.7(2) 


.0010(3; 


.8(1) 


.19(2) 


.3680 


27.8(1) 


0.50(7) 


.22(2) 


14.0(3) 


.0001(3; 


.50(5) 


.22(1) 
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